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Abstract 

I illustrate a general formalism based upon the subtraction method for the 
calculation of next-to-leading order QCD cross sections for any number of jets 
in any type of hard collisions. I discuss the implementation of this formalism in 
a numerical program which generates partonic kinematical configurations with 
an appropriate weight, thus allowing the definition of arbitrary jet algorithms 
and cuts matching the experimental setup at the last step of the computation. 
I present results obtained with computer codes which calculate one-jet and 
two-jet inclusive quantities in photon-hadron and hadron-hadron collisions. 
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1. Introduction 


The study of jet production provides one of the most fundamental tests for the 
predictions of perturbative QCD. The jet rates are rather large, and allow for some 
of the best precision measurements in hadron physics. One-jet and two-jet inclusive 
distributions have been thoroughly measured in the past few years, and an impressive 
amount of data for even more exclusive quantities, like three or more jet cross sec¬ 
tions, has been collected in e + e _ and hadron-hadron collisions. The high transverse 
momenta of the jets set the scale for QCD calculations. The coupling constant is 
therefore small enough to result in a reliable perturbative series. Next-to-leading or¬ 
der predictions for one-jet and two-jet inclusive quantities in hadronic collisions []T]-Q 
have been available for some time. The uncertainty affecting these results is in gen¬ 
eral smaller than the corresponding experimental errors, and a detailed comparison 
between theory and experiments has been carried out. Although some issues need to 
be clarified (like, for example, the tail of the E T distribution measured by CDF ||], 
which appears to be higher than the theory), the overall agreement is quite sat¬ 
isfactory. In the near future, the increased luminosity at HERA will also allow a 
statistically significant study of large transverse momentum phenomena in photon- 
hadron collisions and DIS, thus giving another handle to test next-to-leading order 
QCD predictions QjT . 


The calculation of jet cross sections at next-to-leading order is rather complicated. 
A large number of infrared divergencies is found in the computation of virtual and 
real diagram contributions, due to the large number of colour-interacting, massless 
partons involved in the hard scattering processes. It is then necessary to devise a 
procedure which allows the analytic calculation of the divergent parts and shows their 
cancellation in the sum which defines any infrared-safe physical observable. This task 
has been accomplished in the past by using the slicing H and the subtraction methods. 
In ref. [p] a formalism adopting the subtraction method has been introduced in order 
to calculate one-jet and two-jet inclusive quantities in hadronic collisions. Although 
possible, the extension of this formalism to other type of processes or more exclusive 
observables is not straightforward. For this reason, a fully general formalism based 
on the subtraction method has been proposed in ref. where the results were 
reported in a form motivated by the study of three-jet-like quantities in hadronic 
collisions. The formalism has been subsequently applied to the production of four 
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jets in e + e collisions in ref. 0, and the same calculation techniques have been 
used to study the hadroproduction of jets containing heavy quarks in ref. [|12|] (for 
other applications of the subtraction method, see refs. [[i~3UT7l| ). Afterwards, a general 
formalism has been presented in ref. |l8|] , which uses the subtraction method to cancel 
the infrared divergencies and a new method (called dipole method) to perform the 
analytic treatment of the divergent terms. See also ref. |l9j for another approach. 

The aim of this paper is to illustrate some further improvements of the formalism 


of ref. [ 101, especially relevant for its implementation in a computer code. 1 begin 


by writing the formulae of ref. |T0| in a form such that their generality is apparent, 
and the calculation of any infrared-safe cross section in an arbitrary hard collision 
is straightforward. I then show how to construct a computer code (which I will 
call parton generator) that generates partonic events, which are eventually used to 
plot infrared-safe quantities. In this way, the jet reconstruction algorithm and the 
definition of cuts matching the experimental conditions are inserted at the last step 
of the computation. During the same computer run one can therefore obtain cross 
sections for several, different, jet definitions, as well as predictions for other infrared- 
safe quantities, like for example shape variables. I point out that, in spite of these 
features, such a parton generator is not equivalent to the usual Monte Carlo parton 
shower programs, since it is the result of a fixed-order QCD calculation. Finally, I 
discuss the special case of processes with two or three partons in the final state, and I 
present few results obtained with parton generator coded! written for photon-hadron 
and hadron-hadron collisions, suited for applications to HERA physics. 

The paper is organized as follows: after presenting the formalism in section 2, in 
section 3 I show how to write a parton generator code, and I present numerical results. 
I report my conclusions in section 4. Technical details concerning the formalism 
adopted are collected in appendix A and appendix B. 


2. Formalism 


The goal is to calculate the cross section for some infrared-safe quantity in a 
given hard scattering process. To be specific, I start with the production of IV — 1 
jets in hadronic collisions. According to the factorization theorem in QCD 0 , any 


The codes are available upon request. 
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differential cross section can be written as 


da {HlH2) (K l ,K 2 ) = J2 I dx 1 dx 2 f^ 1 \x 1 )f^ 2 2 \x 2 )dd aia2 (x 1 K 1 ,x 2 K 2 ), (2.1) 


where H\ and H 2 are the incoming hadrons, with momenta K\ and K 2 respectively, 
j s the non-calcnlable but universal distribution function for the parton a, in 
the hadron Hi, and dd aia2 are the (subtracted) short-distance partonic cross sections. 



As shown in ref. |T(J, the cancellation of the infrared divergencies^ arising in the 


intermediate steps of the calculation at next-to-leading order can be fully performed 


at the level of the partonic cross sections. To prove this issue, the fact that initial state 
partons cq and a 2 are quarks and gluons is not crucial. This implies that the same 
proof holds true for quantities like dd ia2 (entering photon-hadron cross sections), dd ea2 
(entering DIS cross sections), as well as for the e + e~ cross section dd e + e ~. For this 
reason, in the following I will only deal with partonic cross sections; it is understood 
that initial state partons will have to be interpreted in a broad sense (that is, they 
can be quarks, gluons, photons and electrons, depending upon the type of physical 
hard scattering process one is interested to study. On the other hand, with final state 
partons I will always mean quarks and gluons). We will get the (N — l)-jet cross 
section in the collision of particles A and B by using the equation 



where (xi, x 2 ) is a suitable luminosity function. Eq. @3) reduces to eq. 0) 
with L^\ X1 ,X 2 ) = ./'if 7 '^ (x-i)fa 2 ^(x 2 ). In the very same way, from eq. (|2.2|) we 
get the photon-hadron cross section if we put 


= '5 7 . 1 '5(1 - xi)/if a) fe), 


(2.3) 


the electron-hadron cross section in the Weizsacker-Williams approximation with 


(4 e) is the Weizsacker-Williams function), the DIS cross section with 

rSd 11 ’ 12 ) = 'Wt 1 - zi)/if 2, 0t2), 


(2.4) 


(2.5) 


3 In this paper, I will never deal with the problem of ultraviolet divergencies. I assume that they 
are renormalized in a proper way. 
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and the e + e cross section with 

L( C.f\ X l, X 2) = Se+ ai S(l - X^Se-aJil - X 2 ). (2.6) 

At the next-to-leading order in QCD, I write the partonic cross sections for the 
production of (N — l)-jets as 


da ai a 2 = da^ 2 + . 

The leading order term dd^\ gets contributions from the processes 


(2.7) 


TV—1 


dl(ki) + a 2 (k 2 ) - ■> ^3(^3) + •'' + a AT+l(^TV+l) ( + e (^e)) j 


( 2 . 8 ) 


where the final state partons a*, i — 3,..., N + 1, are quarks and gluons, the electron 
in the final state is present only in the case of DIS (in which case, cq = e, a 2 = g, q, q). 


1 will denote the processes of eq. (|2.8|) as (N — l)-parton processes, apart from the 
presence of the electron in the final state. We can write 


da^ = _-_ 

“i“2 r_/v — 1)! 


^ 1 ' ) ({az}) V+1 )5Ar_i(i0Ar_i 


(2.9) 


Ddf +1 


In this equation, A4^ N 1 ^({ ci t}i V+ 1 ) is the invariant amplitude for the process of 


eq. ( p.8|) , squared, summed over final state and averaged over initial state colour 
and spin degrees of freedom and multiplied by the flux factor, and d(pN-i is the 
phase-space for iV — 1 massless partons (plus the electron in DIS). In order to in¬ 
clude the contribution of all the partonic processes initiated by cq + a 2 , a sum over 
the flavours g,u,u ,... of final state partons has been performed; the statistical factor 
l/(N — 1)! has therefore to be inserted to avoid double counting. The quantity Sn-i, 
called measurement function, embeds the definition of the momenta of the N — 1 jets 
in terms of the momenta of the N — 1 partons. I will discuss its properties at length 
in the following. 

The next-to-leading order term is 


dd^ = da^ + da ^ + da ^ , 

a^a2 Q'lQ'2 a l a 2 a l a 2 7 


( 2 . 10 ) 
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where 


da[ v l = 


(N — 1)\ In _ xN+ i 


T 1,v \{ai}i +1 )SN-id(j)]sr_-i 

Hdf 


( 2 . 11 ) 


is the contribution of the QCD loop corrections to processes in eq. (|2.8|) , 


da *X = ]v! E M^Ha^+^Sifd^N 


( 2 . 12 ) 


Hdf + 2 


is the contribution of the tree amplitude of IV-parton processes, and da[ c \ is the 


contribution of the initial state collinear counterterms (if needed). Eq. (|2.12|) is anal¬ 
ogous to eq. ( p.9|) . The N — 1 jets have now to be dehned in terms of the IV-body 
partonic kinematics: this is accomplished by the measurement function Sn- 

Although the quantities dehned in eqs. ( j2.11|) and ( p,12|) and the collinear coun- 
terter m s are infrared divergent, their sum in eq. (p,10|) is hnite provided that the 
corresponding observable is infrared safe. As it was shown in refs. j9j.|T(J], this require¬ 
ment can be easily expressed in terms of the measurement functions. Explicitly, the 
conditions (infrared limits) 


lim Sn = <Sn-i , lim Sn = Sn-± , 

feillfc,- 

lim Sjsr = Sn-i , lim Sn = Sn-i , 

fci||fci fc;|l fc 2 


(2.13) 

(2.14) 


with 3 < i < N + 2, 3 < j < N + 2, i ^ j, guarantee that the next-to-leading 
order contribution da^ is hnite. In eqs. (|2.13| ) and ( |2.14|) , Sn-i is constructed 
with the (N — l)-parton kinematics obtained from the A^-parton kinematics in the 
limits indicated. 1 point out that the proof of ref. [|T(]] has been carried out for 
N = A (which corresponds to three-jet production); nevertheless, the fact that N = 
4 has never been explicitly used, and therefore the proof is valid for an arbitrary 
N. Also, the proof of ref. [|TDj] exploits eqs. ( |2.13| ) and ( |2.14| ), which in that case 
follow from the specihc jet definition adopted. Here 1 do not hx a jet definition, and 
therefore eqs. ( [2.13|) and (|2.14|) play the role of conditions which must be fulfilled by 
the jet definition chosen to induce an infrared-safe cross section. Finally, notice that 
eqs. ( |2.13|) and (|2.14|) imply that most of the multiple infrared limits of Sn vanish (by 
definition, in these limits two or more infrared divergencies overlap). This is because 
in these limits the A^-parton kinematics reduces to a configuration where only N — m 
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partons (m > 2) have non-vanishing transverse momentum, and it is not possible to 
define N — 1 jets with less than N — l hard transverse partons. The only non-vanishing 
multiple infrared limits are the soft-collinear ones 


lim __ S N = S N _ i, 

\\kj 


lim ^ S N = S N _ i, 


lim S N = 5iv_i. (2.15) 

k^0,ki\\k 2 


I remark that multiple infrared configurations which are not soft-collinear can con¬ 
tribute to QCD cross sections beyond next-to-leading order. 

In order to prove that eq. ( |2.10| ) is finite, one has to evaluate analytically, with 
some suitable regularization, eq. ( |2.11|) , eq. (|2.12| ) and the collinear counterterms and 
then show that the divergent terms mutually cancel. The structure of the divergencies 
in eq. (|2.11|) naturally arises from the calculation of loop integrals. The case of 


eq. (|2.12|) is more involved: in fact, the divergencies are due to the integration over 
the regions of the phase space where one parton is soft, or two partons are collinear 
(which I will call infrared singular regions). Due to the complexity of the iV-body 
kinematics, one can not perform the analytic integration over the whole phase space. 
The best one can do is to deal with one soft-collinear singularity at a time. To 
achieve this goal in the framework of the subtraction method, in ref. || the iV-body 
matrix elements squared (in that case, N — 3) were decomposed into single-singular 
terms (having, by definition, one soft-collinear singularity at most); each term was 
then integrated over the relevant infrared singular region. Although in principle this 
method can be extended to larger values of N, the amount of algebraic calculations 
and analytic integrations required grows very rapidly, and poses serious difficulties 
already with N — 4. In ref. m a different approach was proposed to overcome 
this problem. The key idea is to use the properties of the measurement functions to 
integrate over the infrared singular regions. In particular, the following decomposition 
can be exploited 

\ 


N +2 


/ 


$n - 


i =3 


+ E S8’»K - kl) 


N +2 


( 1 )/ 


V 


j=3 


(2.16) 


/ 


The terms in the RHS of this equation are defined by their behaviour close to the 
infrared singular regions. In particular 


5,-°^ 7 ^ 0 only if k® —► 0 , /q || k\ , fcj || k 2 


(2.17) 












7- 


5i 1} ^ 0 


y 


only if /q° —> 0, k® —■> 0, /y || k 


(2.18) 


m 


I stress that 2y-°' ) and S^ vanish in the infrared limits not explicitly indicated i 
eqs. ( |2.17| ) and (|2.18|) . Two remarks are in order here. Firstly, eqs. ( [2.171) and (|2.18| ) 
only constrain the infrared limits of and 1 . Therefore, these quantities can be 
redefined up to terms which vanish in these limits. This is the case of the functions 
Sj 1 "' 1 introduced in ref. ||10|| , which have been re-absorbed in the present paper into 
Sj 0 ' 1 and S-j^ (I will show in appendix A that this can be consistently accomplished). 
Secondly, in ref. [10|], eqs. (|2.16| ), (|2.17|) and (|2.18|) have been derived using a given 
jet definition. Nevertheless, as already observed there, it is easy to understand that 
Sj 0 ' 1 and fulfilling these equations can be defined starting from any infrared-safe 


prescription, since they are directly induced by eqs. (|2.13|) and (|2.14|) . I will give 
explicit examples in the following. 

Inserting eq. ([2.16|) into eq. ( p,12| ), and exploiting eqs. (|2.17|) and ( |2.18|) , we see 
that dapJ ao is split into a sum of terms each of which has one soft-collincar singularity 
at most. Therefore, this procedure is equivalent to a single-singular decomposition 
of the matrix elements squared, without requiring any algebraic computation. In the 
end, following ref. [lj]], the result for the next-to-leading order term is arranged as 
the sum of a IV-parton contribution and of a (N — l)-parton contribution 


da {1) = da {1 ’ N) + da^-V. 


(2.19) 


I report in appendix A and appendix B the explicit form of the quantities in the RHS 
of eq. (|2.19| ). Each term contributing to eq. (|2.19|) is finite, and therefore we have an 
operational prescription for the numerical evaluation of an arbitrary jet cross section 
in the framework of the subtraction method. 

3. Numerical calculations 


I now turn to the problem of implementing the formalism of the previous section 
in a computer code. As a benchmark example, I use two-jet production in hadronic 
collisions, and define the jets through the algorithm introduced by Ellis and Soper in 


ref. [21] and formulated in terms of the quantities 


di = kl 


(3.1) 
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= (Vi - Vjf + (<Pi - <Pj) 2 , 

(3.2) 

dij 

= min{k 2 iT ,k 2 jT )^ . 

(3.3) 


The constant D is the jet-resolution parameter. We have ( p d are the jet momenta) 
5 *- 0) = (pj i - %) 5 {pj 2 ~ h) 

cr(J) 

x9(min([di\) - d t )9 (p JlT - 9 ( p JiT - p™ n ) , (3-4) 

= H S (pj i kj ) 6 (pj 2 - h ) 

a(J) 

x9(min([d ij }) - <y 0 (p JlT - ^ <P, 2 t - P^ m ) , (3-5) 


and 5 3 is defined through eq. ( EHP - Here, {i,j,l} = {3,4,5}; H a{J) denotes the sum 
over the permutations of jet labels, with a normalization factor 1/2 inserted, min ( [d ,]) 
(mm([tiij])) is the minimum of the quantities d a , d a p with di (dy) excluded. p™ n and 
p™ n are the minimum observable transverse energies of the two jets (they are fixed, 
input parameters). Finally, when two partons are merged into a jet, eq. ( |3.5|) , the jet 
three-momentum is defined as the sum of the parton three-momenta; other definitions 
would only imply changing the argument of the first 5 in eq. Q3.5I) . When only two 
partons are present in the final state, we have 


Si=J2 s (A. - fe>) (p,, - h) 0 (p v - P’S’ •) 0 (p v - P'S") ■ (3.6) 

a(J) 


It is very easy to check explicitly that eqs. (p.4|), ( fl.5|) and (|3.6| ) fulfill the conditions 
on the measurement functions discussed in the previous section. 

It is now possible to evaluate the quantity 


< H >= £ / HL -„ (da<°l 2 + dai'X) • 

a 1 a 2 


(3.7) 


where L aia2 is the relevant parton luminosity and H is any function of the jet mo¬ 
menta. We may think of H as a product of 9 functions, implementing experimental 
cuts and selecting a bin of a given histogram. < H > will therefore be interpreted as 
the next-to-leading order QCD prediction for the cross section in that bin. The rules 
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for the numerical evaluation of eq. ([h7]) can be read from eqs. ( p.9|) , (|A.3|) , (|A.11|) , 
(|A.16| ) and (|A.26|) . Schematically, one has to perform the following operations: 


1. generate the Bjorken x’s and the partonic kinematics; 

2. evaluate the weight and the jet momenta for this kinematical configuration, 
as specified by eq. O and by the expressions of the partonic cross sections 
in eqs. (|2.9|), ( |A.3| ), ( [A.ll| ), (|A.16| ) and ( |A.26| ). Call an output routine with 
the weight and the jet momenta as entries; 

3. in the output routine, put the weight in the histogram bin selected by the 
jet momenta. 


1 stress that the real computation is actually more complicated than the one 1 outlined 
above, since eqs. (|A.3| ), (|A.11|) and (|A.26| ) require subtractions. 1 will not discuss here 
the Monte Carlo calculation of a subtracted quantity, which is by now a standard 


procedure. The interested reader can find a thorough discussion in ref. JT5|. Following 
the prescriptions implicit in eq. (fT7|), it is possible to write a computer code which 
can calculate < H > for any well-defined quantity H. 1 call this code a jet generator. 
For a given choice of input parameters, it returns event by event the jet momenta (as 
defined by the measurement functions), which are eventually put by the user in some 
histogram bin. 


By using the formulae collected in appendix A and B, the construction of a jet 
generator is therefore fairly straightforward, and allows the computation of the cross 
section for the production of any number of jets in any hard collision in the framework 
of the subtraction method, without requiring algebraic manipulations of the partonic 
transition amplitudes, as in ref. ||. The measurement functions, embedding the jet 
definition, can be written by the user in his own computer routine. 


The main drawback of such a jet generator is the following: the jet definition is 
used, through the measurement functions, to disentangle the singularities appearing 
in the real contribution. Therefore, to get the jet cross section for several, different, jet 
definitions, one has to perform several computer runs. Furthermore, a jet generator 
outputs only jet momenta, which implies that it is not possible to calculate non-jet- 
like infrared-safe observables, as for example shape variables. 1 will now show that 
it is not difficult to overcome these problems using the formalism of ref. [ID]. To be 
specific, I start by considering again two-jet production is hadronic collisions. With 
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three partons in the final state, I introduce the following quantities, defined in terms 
of the partonic momenta 


pf 1 = e ,< 0) e(k 3T + k 4T + h T - e™"), 


(3.8) 


'pW 


= 0 !? 0(hr + k lT + k 5T - E"“ n ), 


p 3 = £ 


»=3 


\ 


3=3 

I jW 


(3.9) 

(3.10) 


where 0® and 0^ are suitable products of 9 functions. With two partons in the 
final state, I define 

V 2 = l -9{k 3T + h T - E™ in ). (3.11) 

1 require that, close to the infrared singular regions, the quantities defined in eqs. (|3.8|) 
and ([3T9| ) have the following properties 

7 ^ 0 only if k® —> 0, k t || k\, k t || k 2 , (3-12) 

V\f 7 ^ 0 only if /q° —> 0, > 0, k t || kj (3.13) 


(which means that they vanish in the infrared limits not explicitly indicated), and 
that they are such that (a relabeling of the partons is understood, if necessary) 


X 

lim Vz 

fci||fcl ki\\k 2 

Furthermore, the following equation must be fulfilled 


V2, 

lim P 3 = V 2 , 

ki\\kj 

(3.14) 

V2, 

lim P 3 = V 2 . 

(3.15) 


/ 


E 

i=3 


ei 0) + £e« »(*£.-*; 




4 ) 


V 


3=3 


= 1 . 


(3.16) 


/ 


If we choose the free parameter E™ m such that E™ in < 2 min (p™ m ,p™ n ), e q. ( |3,16| ) 
implies that 

^3 = S%P 3 , S 2 = S 2 V 2 , (3-17) 


where the S functions were defined at the beginning of this section (notice that it is 
always possible to choose E™ in without any reference to a specific jet definition. For 
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example, E™ m can be less than twice the minimum transverse momentum observable 
by the detector). I now define quantities analogous to the partonic cross sections 
appearing in eqs. 


and (|2.19|) (see also appendix A), by formally substituting the 
S functions with the V functions. 1 adopt the following notation 


dirf^V 2 ) 

- da[°l 2 \s 2 ^v 2 , 

(3.18) 


= da {l 

a l 

a 2 1 S 2 -^V 2 1 

(3.19) 

d^d>(P, <0> ) 

= do[ “•/> | sfwr . 

(3.20) 


= d<T i ia J] ,> 

(3.21) 


5 



ksp 3 ) 

= E 

d^(P! 0> ) + Zd^(P!j > ) 

(3.22) 


3 =3 

V jA* / 


The quantity 


[i/5] = £ / L, m ( HS 2 d*<°>JP 2 ) + HS, + HS 3 da™(P 3 )) (3.23) 


can be numerically evaluated, being finite. This can be easily understood by observ¬ 
ing that eqs. ( |3.18|) - (|3.22|) are finite thanks to eqs. (|3.12|) - (|3.15|) ; the proof is identical 
to the proof of ref. |T0 |, which showed that eqs. TO , ( TO , ( CT ) and CT) guar¬ 
antee the finiteness of the next-to-leading order cross section. Furthermore, thanks 
to eq. ( |3.17| ), we have 

<H>= [HS ], (3.24) 

where < H > was evaluated in eq. ( |3.7p . Therefore, every quantity which can be 
calculated with eq. (|3.7|) can also be calculated with eq. ( [3.231) . Nevertheless, eq. (|3.23| ) 
is much more flexible. In fact, its numerical evaluation requires the following steps 


1. generate the Bjorken x's and the partonic kinematics; 

2. evaluate the weight for this kinematical configuration, as specified by 
eq. ( p.23| ) and by the expressions of the partonic cross sections in eqs. ( |3.18|) - 
(|3.22|) . Call an output routine with the weight and the parton momenta 
as entries; 
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3. in the output routine, define the jet momenta as specified by the jet¬ 
finding algorithm embedded in the S functions, and put the weight in the 
histogram bin selected by the jet momenta. 


Therefore, the code implementing eq. (|3.23|) returns event by event the parton mo¬ 
menta, which are eventually manipulated by the user. I call such a code a parton 
generator. The main difference between a jet generator and a parton generator can 


be read from eqs. (|3.7|) and (|3.23|) ; while in eq. (|3.7| ) the measurement functions S 
enter the partonic cross sections and render them finite, in eq. (|3.23| ) they have the 
same role of the function H. In eq. (|3.23|) the partonic cross sections are finite tha nk s 
to the V functions, which can be chosen once and forever without any reference to a 
specific jet definition. For this reason, such a parton generator is able to plot, on an 
event-by-event basis, any infrared-safe quantity defined with two or three partons in 
the final state. Namely, in a single run it can produce one-jet and two-jet observables, 
with the jets defined by several algorithms, as well as various shape variables. Notice 
that eqs. ( |3.11| ) and ( p.lG| ) guarantee that the events are generated in the whole phase 
space. Some of them are eventually rejected, namely those having a kinematics which 
is not fulfilling 


0{k 3T + &4t T ^5t — E™ in ), 
0(k 3T + k AT - E™ n ). 


(3.25) 

(3.26) 


Nevertheless, E™ in can always be chosen in such a way that eqs. ( 3.25 ) and ( 3.261 ) 
are less stringent than the physical cuts which are required to define any infrared- 


safe quantity. Technically, eqs. ( 3.25 ) and ( 3.26 ) are inserted to avoid multiple non 
soft-collinear infrared singularities, which do not contribute to the cross section at 
next-to-leading order. In the soft-collinear limits we have, analogously to eq. (|2. 15 ), 


lim V 3 = V 2 , lim V 3 = V 2 

k^oSSj k°^o,kSi 


lim V 3 = V 2 . (3.27) 

fc0^0,fei||fe2 


It should be clear that the method described so far can be used, without any 
modification, to write a parton generator for photon-hadron or e + e _ collisions. In the 
latter case, we can also set E™ in = 0, since multiple soft singularities are forbidden by 
energy conservation and configurations with partons collinear to the incoming leptons 
are not singular. 
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The situation is slightly more complicated with more than three particles in the 
final state in DIS (two or more jet production), and for three or more jet production in 
hadron-hadron, photon-hadron and e + e~ collisions. Indeed, it is easy to understand 
that with a four-body (or more) kinematics, the analogous of eqs. ( [3.25D and (|3.26|) are 
not enough to avoid multiple non soft-collinear infrared singularities. Eq. ( |3.25|) and 
eq. (|3.26|) have to be substituted by suitable products of 6 functions which vanish in 
the multiple non soft-collinear infrared regions. In the end, one should get quantities 
Vn and Vn- i, defined analogously to Vs and V2 in eqs. ( |3.10|) and (|3.11|) , which fulfill 
eqs. (|3.12| )- (|3.16|) and (|3.27|) (with the formal substitutions 2 —> N — 1, 3 —> N). The 
products of 6 functions which substitute eqs. ( |3.25|) and (|3.26| ) must be chosen in such 
a way that 

Sn = SnVn , = Sn-iVn-i • (3.28) 


As already observed before, I point out that, in the case of hadron-hadron, photon- 
hadron and e + e _ collisions with two or three partons in the final state, with a single 
choice of the V functions it is possible to fulfill eq. (|3.28|) for any measurement func¬ 
tions S. The same is not true for the other cases. In practice, if one chooses the 
V functions in such a way that they only vanish extremely close to the multiple non 
soft-collinear infrared regions, then eq. (|3.28| ) holds for all the physically meaningful 
choices of the measurement functions S. We therefore get again a parton generator 
which is able to plot, event by event, (N — 1)-jet inclusive quantities with several dif¬ 
ferent jet definitions, and shape variables which get contributions at next-to-leading 
order from (N — 1)- and Wparton configurations. 


In order to write a code for a parton generator, a definite choice for the V functions 
has to be made. Although these functions can be rather freely chosen, the only 
constraints being eqs. (|3.12| )- (|3.16|) and ( |3.27|) , their properties are motivated by the 
properties of the measurement functions of a jet algorithm. Therefore, it is quite 
natural to use a jet algorithm to construct the V functions. Notice that, by definition, 
this jet algorithm has nothing to do with the jet algorithm(s) eventually used to obtain 
predictions for physical observables. 


I now consider the case of photon-hadron and hadron-hadron collisions with two 
or three partons in the final state. Using the prescription of ref. |U|, we get, as can 

































be also directly seen from eqs. ( p.4|) and ( fh5p , 


(3.29) 

(3.30) 


@,| 0) = 6{inin([di}) — di), 
= 9(min([dij}) — dij), 


where the relevant quantities have been introduced in eqs. (H)-®' 
fulfill eq. ( |3.16|) , we must have D < 2it/3. If we adopt the prescription 
algorithm |22|], we get 


In order to 
of the cone 


0- O) = 9(Rij - gij)9(Ru - gn)9(Rji - <| j i)^(min(fc iT , k iT ) - k iT ), (3.31) 
@p ) = 9(gij-R ij )9(Rii-gu)9(R j i~ g j i)9(ki T -min(k iT ,k jT )), (3.32) 


where 


9ij 


Eq. (p.l6|) requires that R < tt/3. 


k iT + k 


•jT 


max(ki Tl kj T ) 


R. 


(3.33) 


I wrote a fortran code for a parton generator in photon-hadron collisions, and a 
fortran code for a parton generator in hadron-hadron collisions. The latter can be 
used also in the case of hadronic photon-hadron collisions, that is for photon-hadron 
interactions where the photon fluctuates into a hadronic state before undergoing a 
hard collision. Therefore, the codes are suitable for applications to HERA physics. In 
order to test the codes, I produced single-inclusive jet and two-jet observables. The 
results obtained with the photon-hadron (hadron-hadron) code have been compared 
with the results of ref. || (refs. (6j.[|). In both cases, I found nice agreement. 

In figures |l] and |2] I show the transverse energy of the single-inclusive jet, the 
azimuthal distance and the invariant mass of the pair of the two hardest jets, and the 
transverse thrust distributions for ep collisions in the HERA energy range, E CM = 
300 GeV, as predicted by the aforementioned parton generators. The Weizsacker- 
Williams approximation has been adopted (the form of ref. [^3j, which includes non- 
logarithmic singular terms, has been used with Ql vw = 4 GeV 2 and 0.2 < y < 0.8), 
and therefore we are dealing with a photoproduction process. As such, both the 
pointlike photon cross section (whose contribution, obtained with the photon-hadron 
code, is shown in figure |l]) and the hadronic photon cross section (figure obtained 
with the hadron-hadron code) are sizeable. I point out that the curves presented in 
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Figure 1: Jet observables and transverse thrust in ep collisions 

(Weizsacker- Williams approximation) at HERA. Pointlike photon only. 
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Figure 2: Jet observables and transverse thrust in ep collisions 

(Weizsacker- Williams approximation) at HERA. Hadronic photon only. 


the figures are not physical quantities, if taken separately, since only their sum is 
measurable. Nevertheless, here I am only interested in discussing the properties of 
parton generator codes. Phenomenological results for HERA physics will be reported 


in a forthcoming publication |2J|. The four plots appearing in the same figure have 
been obtained with a single computer run. The jets have been defined using the 
cone algorithm |22|] with R = 1 (solid histograms), and the algorithm of ref. [21 


with D = 1 (dashed histograms). In the same run, I also computed the jet cross 
sections with the cone algorithm and R = 0.7. I do not show the results, since 
they are extremely close to the dashed histograms. I used the set MRSA' for the 
parton densities in the proton, and the set GRV-HO for the parton densities in the 
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photon. The value of A QCD was fixed at the value suggested by the MRSA' set, 
A 5 = 152 MeV. Factorization and renormalization scales have been taken equal to 
half of the total transverse energy of the event. The figures have been produced using 
eqs. ( [3.31| ) and (|3.32|) , with R = 1 , to define the V functions; 1 verified that completely 


equivalent results can be obtained by using eqs. (|3.29|) and (|3.30| ) with several values 
of D. This gives a consistency check on the method used to construct the parton 
generator codes. Also notice that the curves corresponding to the two different jet 
definitions adopted have comparable quality. Naively, one may think that the curves 
relevant for the jet definition which matches the prescription for the V functions could 
converge somewhat faster. This is not true; the convergence properties are almost 
completely dictated by the behaviour of the V functions and of the jet definitions close 
to the infrared singular regions, where they are strongly constrained by the infrared 
safeness conditions, which are definition-independent and must hold in any case. 

Finally, 1 would like to comment on the numerical calculations performed with 
the subtraction method. The main drawback of the method is that in the subtracted 
integrals the cancellation takes place between terms which have different kinematics 
(see for example eq. (|A.5|) ). For this reason, in ref. [1_0;] the arbitrary parameters £ cu4 , 
8 U 5 0 have been introduced in the calculation (see also refs. |L5|-[n|]). To be specific, 
consider eq. (|A.5|) ; we see that the counterevent (that is, the integrand function 
calculated on the pole, /(0)) is subtracted only if £, < £ cut . Therefore, with a suitable 
choice of £ cut , the counterevent is subtracted only if the event is close to the soft 
limit £j = 0 . In this way, the subtraction is performed only in those cases when 
the kinematics of the event and of the counterevent are quite close to each other, 
thus resulting in an improved numerical stability of the result. As a bonus, we 
save computing time, since in most of the cases the calculation of the counterevent 
is not necessary. It is very important to realize that the parameters f cut , 8 I} S a 
have nothing to do with the non-physical parameters which must be introduced in 
the slicing method. In the present case, the parameters do not have to be small, 
since no approximation was performed in the intermediate steps of the calculation. 
Therefore, there is no need to prove that the physical cross section is independent of 
icut-, 5 0 , since this is true by construction. We can exploit this fact as a powerful 
test of the numerical implementation of the formalism, by verifying that the cross 
section is a constant with respect to the choice of f cut , 5 T , 5 a . By choosing ( f cut , 8 r , S Q ) 
close to ( 0 , 0 , 0 ), we also see that the slicing method can be obtained as a limit of 
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the subtraction method (for this to be rigorously true, the kinematics of the events 
close to the infrared limits must be set equal to the kinematics of the corresponding 
counterevents. For example, this amounts to set /(£,) = /(0) for & < icut ineq.flAJ)). 
In this case, as it is apparent from eqs. ( |A.3|) , ( |A.11|) and (|A.16|) , cancellations between 
large numbers take place in the sum which defines the physical cross section. 


Conclusions 


In this paper I have studied the definition of infrared-safe cross sections at next- 
to-leading order in QCD. The formalism of ref. [|l0j, which is based on the subtraction 
method, has been written in a form which clearly shows its generality and universality. 
The formulae can be applied to any hard scattering process, with an arbitrary number 
of final state partons. I then turned to the problem of the numerical computation of 
infrared-safe cross sections. It has been shown how to write a computer code which, 
event by event, outputs the momenta of the final state partons. The momenta are 
eventually used in an analysis routine to define the physical observables. The special 
case of infrared-safe quantities defined with two or three partons in the final state has 
been studied. Two codes, one dealing with photon-hadron collisions and one dealing 
with hadron-hadron collisions, have been written. Sample numerical results have been 
presented for the case of ep collisions in the Weizsacker-Williams approximation. This 
is relevant for applications to HERA physics, where both the pointlike photon and the 
hadronic photon give sizeable contributions to the physical cross section. I presented 
predictions for one-jet and two-jet observables (as an example, I defined the jets using 
two different prescriptions), and for transverse thrust. 
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APPENDIX A: N-parton contribution 


The IV-parton contribution, first term in the RHS of eq. (|2.19|) , is written as the 
sum of finite terms, whose form is induced by the decomposition of the measurement 
function of eq. (|2.16| ) 


N+2 


( 


da = y 

“1»2 / 


1=3 


N+2 


\ 


da {inJ) + y da (out ’ f) 

a l a 2 Z—✓ a l a 2> l 3 


\ 


3 =3 

j¥=i 


(A.l) 


The hrst term in the RHS of this equation is proportional to <S,- 0 , the second one to 
Sjj. By recalling the properties of the measurement functions close to the infrared 
limits, it is easy to choose a set of variables suitable to perform the calculation with 
the subtraction method. In the partonic center-of-mass frame, we parametrize the 
momentum of the final state parton i as 


ki 



(A.2) 


where a/A is the partonic center-of-mass energy. We get |10 


da {inJ) = 

a-^a 2 ,1 


x 




+ 


<5/ 




(1 - y Af (iV) ({cq}f +2 ) s\ 0) dcjxiZidyidipi . 

V {««> 3 r+a ) 


(A.3) 


The measure d(f> is implicitly defined by writing the phase space for N partons (plus 
electron in DIS) in 4 dimensions as 


d(j) N 




2 

iidiidyidipi. 


(A.4) 


The statistical factor 1/N\ is due to the sum over the flavours of final state partons 
(see eq. (|2.12|) ). The damping factor (1 — yf )tf and d^ 0 " 1 guarantee that the integrand 
is finite everywhere in the phase space. We can therefore use 



, .fit;)-f (mu, -6) 


/ o 


(A.5) 
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(A.6) 




g{y%) -g(±i)0(±yi 
It Vi 


1 T 5 : ) 


which hold for any smooth functions / and g. These prescriptions are almost identical 
to the usual T prescription, except for the fact that the value of the integrand function 
at the pole is subtracted only if the integration variable satisfies the condition imposed 
by the 6 functions. The parameters 


0 < icut < 1 , 0 < 8j < 2 


(A.7) 


can be arbitrarily chosen in the indicated range. Finally, notice that if vanishes 
in the soft (£, —> 0) and collinear (r/j —> ±1) limits, eq. (|A.3| ) becomes 




(A.8) 


{«df +2 


This expressions is identical to the one induced by the functions S[' n in ref. [[H] 


and therefore proves that the contribution of S{ m can be safely re-absorbed into the 


contribution of . 

In order to give the explicit expression for the second term in the RHS of eq. (|A.1| ), 
we rewrite eq. (|A.2|) as 




h = —& (l, k) , k = piR , ^ = (O, l) , 


(A.9) 


where A is a suitable matrix, and parametrize the momentum of parton j as 

Vs. 


kj = 

J 2 


Tj f 1 ,T) . k,=P , R . = 


(A. 10) 


Therefore, in the limit yj —> 1 the partons i and j become collinear to each other. 
We get |L0 


da {out ’ f) = — ( -) ( E_ 
N\ \&J C V 1 - Vj, 


x Sf^(k) T -kl) 


( 


(!-%■)*& E 


} 3 " +2 


'VS' 


2(2tt) 


dcftd^id^jdyidyjdipidpj , (A. 11) 
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where 


[ dyjg{yj )[—^—) 


f 1 d g(yj) - gO-Mvj -1 + £>) 

•/-I J 1 - Vj 


(A-12) 


and 0 < 5 a < 2. The phase-space for A partons (plus electron in DIS) in 4 dimensions 
has been written in the following way, thus implicitly defining dp 


d(j)N — dtp 


i (Vs' 

2(2vr) 3 l 2 / 


2 \ 2 


d£i d^jdyi dyjdipid^j 


(A.13) 


If sV vanishes in the soft (£* —> 0) and collinear (■yj —> 1) limits, we get 


KZ’Ihpr E - *£•)#»■ (A.i4) 


Rdf +2 


Like eq. (|A.8|) , this shows that the contribution of <S/ m can also be re-absorbed into 


the contribution of <S,'j ■*. 


1.1. (N-l)-parton contribution 


We now turn to the second term in the RHS of eq. ( j2.19| ). We write it in the following 
way 

da ( iA-i) = da (1 ’ N - lv) + da (1 ’ N ~ lr) , (A.15) 

a l a 2 a l_ a 2 a l a 2 ' \ / 

where the first term reads 


K’t 1 "' = C E S(Mf +1 )*r (0) ({a i }p 1 ) 


2tt 


Rdf +1 


- i Ar + 1 

+ E A"’({“i}f +1 ) 


( 0 ), 


2tt 2 


*»j= 1 

At 


Rdf^ 
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27 t (A 1 )! , , w+1 


Rdf 


(A.16) 
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Here 


J 


(reg) 


U 


87T 2 


1 log 2 + log ^ log k ' 3 ' kl 

2 8 Q 2 8 g 2 K 2 Ej Ej 


2 kj ■ ki 


2k.j ■ ki ' 


/ krh_ 

\2E j E i/ 

ki ■ h 


log 2 ' — log 4- 2 ' log 2 1 — 2 log 2 2 


EjEj 


EjEi 


2EjEi 


and 

Q(Mf +1 ) 


N +1 


= E 

J=3 


7 '(«j) - log 7777 (7 K) - 2 C'(a j ) log ■ 2jE; ' ^ 


2Q 2 V' V “ J/ ~~^'~°LutVs) 


+ 2C{dj) log' 


2 2 -% 




log 2 £cut - 27 (aj) log 


2 -^j 


(A.17) 


- log ^ 7 ( a i) + 2C(ai) log C, cut + 7 (a 2 ) + 2C(a 2 ) log j . (A. 18) 

Here £) is the energy of parton i in the partonic center-of-mass frame, /i is the 
factorization scale, and Q is the Ellis-Sexton scale (see below). The expressions for 
the colour factors C(a), 7 (a) and 7 '(a) are 


C{g)=C A , C(q) = C F , 

.. 11C A — AT F Nf 
7 (g) = -g- £ , 


7(9) = 2 Cf ’ 


„ , 67^ 27t 2 23 Ar ,, , 13 

7 (fl 1 ) ~ ”g _( ^' 4 ~3~^ a “cf^E ’ 7 (?) — ~2^ F 


2vr 2 


a. 


(A.19) 
(A. 20 ) 

(A. 21 ) 


Notice that eq. (|A.18l) depends upon the flavours of the initial state partons. In the 
case of photon-hadron collisions (DIS), we have cq = 7 (ai = e), and in the case of 
e + e _ collisions we have cq = e + , a 2 = e~. Eq. (|A.18|) still holds, and we just define 


C(a) = 0 , 7 (a) = 0 if a — 7 , e . 


(A. 22 ) 


We also defined 

dcr(°)({a i }f+ I ) = —l- I y i X<' v - 1 >({a,}f'+ 1 )S N _ 1 # N _ 1 , (A.23) 
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which is identical to eq. (|2.9|) except for the fact that the sum over the flavour of final 
state partons is not performed, and 




(N — 1)!' 


1 \{ a i}i +1 )^N-id4>N-i- (A.24) 


The functions are usually denoted as colour-linked Born squared amplitudes. 

They enter the expression of the virtual corrections to the (N — l)-parton processes, 
eq. (|2.11|) . We adopt the following form p5| .j4| 


M (n- i,v) = as ( 4 tt) £ ffi 2 ' 

27T T(1 — e) \Q 2 y 


i N+l 1 N+l \ 

- ^ C((ii) + -J2 7(o<) ].M (iV-1) 


i— 1 


2=1 


+ — Y log 

2e ^ g Q 2 8tt 2 lJ ^ 


i,j =1 

*7 


(A.25) 


This equation also defines the finite part of the virtual correction, Ai^ s 1,1; \ which 
enters eq. (|A.16|) . 

The second term in the RHS of eq. (|A.15|) is 


dd (1 ’ N ~ lr) 
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dCL2 


x dcr (0) (ai,d, {ai}^ +1 ;k u (l - t,)k 2 ) d£, (A.26) 

(«df +1 


where P^(z, 0) + eP^f(z, 0) + 0(e 2 ) is the Altarelli-Parisi kernel for z < 1 in 4 — 2e 
dimensions, the sum J2 d runs over g,u,u ,..., K a j, define the finite part of the initial 
state collinear subtraction (in the MS scheme they are equal to zero), C a b = 1 when 
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b is a quark or a gluon and, analogously to eq. (|A.5|) , 


In the case of photon-hadron collisions, we define 

P dj (z) = 5 dq ^P qg (z), C dl = S dq e 2 q ^ . (A.28) 

The functions K dr , depend upon the scheme adopted for the partonic densities in the 
photon, entering the hadronic photon-hadron cross section (for details on this issue, 
see for example ref. |^6[). Finally, P de = K de = 0, C de = 1. This formal statement 
corresponds to the fact that no collinear singularity is associated with an incoming 
lepton. 


/(O-/(o)0(6*t-O 


log£ 


(A.27) 








